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Spatially unassociated galaxies contribute significantly to the 
blended submillimetre galaxy population: predictions for follow-up 
observations of ALMA sources 
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ABSTRACT 

There is anecdotal evidence that spatially and physically unassociated galaxies blended into 
a single submillimetre (submm) source contribute to the submm galaxy (SMG) population. 
However, the significance of this subpopulation has neither been observationally constrained 
nor theoretically predicted. This work is the first to theoretically predict the contribution of 
spatially unassociated components to the SMG population. We generate mock SMG cata- 
logues using lightcones derived from the Bolshoi cosmological simulation; to assign submm 
flux densities to the mock galaxies, we use a fitting function previously derived from the re- 
sults of dust radiative transfer performed on hydrodynamic al simulations of isolated disc and 
merging galaxies. We then calculate submm number counts for different beam sizes and with- 
out blending. Our model suggests that there are a sufficient number of blended SMGs to ac- 
count for the observed number counts of submm sources with 850-|im flux density 5*850 ^12 
mJy. Furthermore, we predict that ^30 — 80 per cent of blended SMGs have at least one spa- 
tially unassociated component. For a 15-arcsec beam, blends of > 2 galaxies in which at least 
one component is spatially unassociated dominate the blended sources with total 5*850 ^ 3 
mJy. The distribution of the redshift separations amongst the components is strongly bimodal. 
The typical redshift separation of spatially unassociated blended sources is ~ 1. Our predic- 
tions for the contributions of spatially unassociated components and the distribution of redshift 
separations are not testable with currently available data, but they will be easily tested once 
sufficiently accurate redshifts for the individual subcomponents (resolved by, e.g., ALMA) of 
a sufficient number of single-dish-detected blended SMGs are available. 

Key words: galaxies: abundances - galaxies: high-redshift - galaxies: luminosity function, 
mass function - cosmology: theory - cosmology: large-scale structure of Universe - submil- 
limeter: galaxies. 



1 INTRODUCTION 



Subm illimetre (submm) galaxie^j (SMGs; see iBlain et alj 
120021 for a review) are some of the most luminous. 



* E-mail: christopher.hayward@h-its.org 

f CITA National Fellow 

1 Throughout this work, we use the acronym SMG in the conventional 

manner - to denote single-dish-detected submm sources - although a given 

SMG may be composed of multiple distinct galaxies. This issue is discussed 

in Section|4.2| 



rapidly star-forming galaxies (wi th star forma t ion r at es of 
~ 10 2 - 10^ Mb yr _1 ; e.g., iMagnelli et all |2010| . 120121 : 
IChapman et ail 120100 in the Universe. Although simulations 
have suggested that merger-induced starbur sts can produce 



the s u bmm fluxes cha r acteris ti c of S MGs ( Chakrabarti et al. 



2011). 



I2008J : Narayanan et"ail 120091 |2010| : lHavward et al 

explaining the observed number density of such prodigious 
star-formers has posed a challenge for semi-a nalytic mod- 
els a nd hydrodynamical simulations ( e .g.. iGranato et al 
20001 : iBaughetalJ 120051; iFontanot et all l2007f: iDave et al 



20 id : iNaravanan. Bothwell. & Pavel I20I2T Narayanan & Dave 
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20121: iNiemi et all |2012|; IShimizu. Yoshida. & Okamotd I2012J : 



Somerville et al.l 12012c lHavward et al.l 120131) . which has ca used 
some to consid er them a challenge fo r ACDM <|Primackll2012l) . 

Recently, lHavward et ail 1201 ll |2012|. l2013h suggested that 
treatment of the effects of blending caused by the large beam sizes 
(FWHM - 15 arcsec, or ~ 130 kpc at z ~ 2 - 3) of the 
single-dish telescopes used to identify SMGs can help alleviate the 
discrepancy between the SMG number counts predicted by mod- 
els and those observed. Specifically, during the pre-coalescence 
phase of major mergers, b oth progenitor discs can be blended 
into a single submm source: lHavward et al.l ( I2013L hereafter, H13) 
suggest that such 'galaxy-pair' SMGs account for a significant 
fraction (> 30 per ce nt) of the SMG population . Interferomet- 
ric observations (e.g., Tacconi et al. 2006, 2008: Bothwell et al. 



201CtlRiechers et alfeo 11 alibi: lEngel et alj|20ld:lwang et alll201ll : 



Barger et al.l 120121 : iKarim et al.l 120121 : ISmolcic et all 120121) have 
provided many examples of single-dish SMGs that are resolved 
into multiple components when observed at an order-of-magnitude 
better resolution, and in some cases, the redshifts and kinematics 
suggest that the individual sources are widely separated (projected 
separation > 10 kpc) discs in the process of merging. Recent inter- 
ferometric continuum imaging surveys suggest that of order half of 
single-di sh-detected SMGs are blends of two or m ore distinct com- 
ponents iKarim et alj|2012l : ISmolcic et alj|2012l) . Interestingly, in 
those surveys, the brightest sources (those with 870- pm flux den- 
sity 5*870 > 12 mjy) are almost all blended sources. This cutoff 
suggests an upper l imit on the star formation rate (SFR) of SMGs 
JKarim et alj|2012l ; lHavwardl2013l) . which may be a consequence 
of feedback from star formation and/or active galactic nuclei but 
may also be simply a consequence of limited gas supply. It also 
suggests that to reproduce the brightest SMGs, models must ac- 
count for the effects of blending. In this work, one question that 
we address is whether we can reproduce the number counts of the 
brightest SMGs via blended sources alone. 

Early-stage mergers are only one type of blended submm 
source; it is also possible that spatially (and thus physically) 
unassociated projected multiples contribute to the SMG popula- 
tion. There is alre ady anecdotal evidence for this subpopulation: 
IWang et alj u20 1 lh present two examples of SMGs that are resolved 
into multiple components by the Submillimeter Array (SMA) and 
for which the resolved components are located at significantly dif- 
ferent redshifts (in one case, the difference between the two com- 
ponents is Az = 1.015, and in the other, the three components 
have redshifts of 2.914, 3.157, and 3.46). Thus, the individual com- 
ponents are not widely separated galaxies that will merge but rather 
chance projections of completely unrelat ed galaxies. Some of the 
sources observed bv lSmolcic et al.l a2012b may also be examples of 
this subpopulation. 

The relative contribution of such projected multiples to the 
SMG population is currently unconstrained. However, this situation 
is likely to change soon, once redshifts for a significant number of 
the sources from the Atacama Large Millimeter/submill imeter Ar- 
ray (ALMA) follow-up observations JKarim et alj|2012l ; Hodge et 
al, submitted) of t he Large APEX Bolometer Camera (LABOCA; 
ISiringo et a l. 2009) Exte nded Chandra De ep Field South Submil- 
limetre Survey (LESS: IWeiB et al.l 120091) sources are available. 
Thus, a prediction for the relative contributions of spatially asso- 
ciated and unassociated components is particularly timely. 

In this work, we make such a prediction. To do so, we use halo 
catalogues derived from a cosmological simulation. Stellar masses 
and SFRs are assigned using observationally constrained stellar 
mass- and SFR-halo mass relations. For each galaxy, we assign an 



850-|i.m flux density Ssso using a fitting function, which was previ- 
ously derived from the results of dust radiative transfer performed 
on hydrodynamical simulations of isolated disc and merging galax- 
ies, that gives Ssso as a function of SFR and dust mass. Then, we 
search the mock catalogues for blended submm sources and calcu- 
late the total submm flux density for each blended source. Using a 
cut based on the redshift separations of the individual components 
of the blended sources, we separate the blended sources into 'spa- 
tially associated' and 'spatially unassociated' subpopulations and 
calculate the relative contributions to the submm number counts as 
a function of single-dish submm flux density. 

The remainder of this work is organised as follows: in Section 
2, we discuss the details of our method. In Section 3, we present our 
predictions for the cumulative submm number counts for different 
beam sizes, the relative contributions of spatially associated and 
unassociated sources, and the distribution of redshift separations 
amongst the individual components of blended sources. In Section 
4, we discuss some implications of our results. Section 5 presents 
our conclusions and plans for future work. 



2 METHODS 

2.1 Creating mock galaxy catalogues 

The method used in this work differs significantly from that of H 1 3 
because the method used in H13 cannot treat spatially unassociated 
projected multiples. Here, we use mock lightcones derived from 
a cosmological simulation to provide information about the large- 
scale clustering of dark matter haloes. We start with halo catalogues 
from the Bolshoi simulation iKfypin. Truiillo-Gomez. & Primackl 
l201ll : lBehroozi et al.l2013al lbl). which is a collisionless dark matter 
simulation performed using the Adap tive Refinement Tree (ART; 
iKravtsov. Klypin. & Khokhlo\lll997l ) code. The Bolshoi simula- 
tion volume has side length 250 h^ 1 Mpc and contains 2048 3 par- 
ticles, which yields a mass resolution of 1.35 x 10 s h~ 1 Mq. The 
force resolution is 1 h" 1 kpc (physical). The cosmological param- 
eters used for Bolshoi are Q m = 0.27, Qa — 0.73, h — 0.7, 
erg = 0.82, and n s — 0.95. Th ese parameters are co nsistent with 
the WMAP7+BAO+Ho results JKomatsu et"ai]|201ll) . The values 
for us and n s are also consistent with the recent Planck results 
jPlanck Collaboration et al.l2013al lbl). The value of Q, m (Q,a and h) 
used is 2.6a greater (2.2a less) than that from Planck. Neverthe- 
less, the uncertainty in the cosmological parameters should affect 
our results less than the other uncertainties discussed below. 

Starting at eight random locations within the simulation vol- 
ume, we select haloes along a randomly oriented sightline with an 
84' x 84' (1.96 deg 2 ) field of view from z = 0.5 to z = 8; this 
redshift range contains the bulk of the SMG population, both in the 
real U niverse (e.g. JSmo lcic et al. 2012; Yun et al. 2012; WeiB et al. 
120131) and in our model (for reasons discussed below). Cosmolog- 
ical redshifts are calculated including the effects of halo peculiar 
velocities. Stellar m asses are assigned based on the s tellar mass — 
halo mass relation of lBehroozi. Wechsler. & Conrovl 02012b . which 
accounts for redshift-dependent scatter in stellar mass at fixed halo 
mass and the effects of systematic observational biases in stellar 
mass reco very. SFRs are assigne d based on the SFR — halo mass 
relation of lBehroozi et all ( 120121) : this includes redshift-dependent 
scatter in SFRs at fixed halo mass and the increasingly tight cor- 
relation between stellar mass and SFR at fixed halo mass at high 
redshifts. To approximate environmental effects near massive clus- 
ters, satellite SFRs are reduced by a factor equal to their current 
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Figure 1. Demonstration of the difference between spatially associated and spatially unassociated galaxy pairs. Open circles denote the locations of mock 
galaxies, and dashed lines demarcate cubic regions with side length 10 h^ 1 Mpc (comoving). The focus is on galaxies in the central cube and their neighbours 
(which may be located in the neighbouring cubic regions). Left: X-Y plane. Red solid symbols denote galaxies that are associated in three dimensions, 
whereas blue solid symbols denote galaxies that appear to be associated in projection but are actually unassociated in three dimensions. An observer without 
access to accurate redshifts would be unable to distinguish between these two cases. Right: X-Z plane. This rotated perspective verifies that the galaxy pairs 
denoted with red symbols are associated, whereas those denoted with blue symbols are certainly not. The lengths of the cyan lines connecting the unassociated 
galaxy pairs are substantially larger than their apparent separations in the projected X — Y plane. 



subhalo mass divided by the peak mass in their subhalo's mass ac- 
cretion history. This is only a toy model for satellite quenching, of 
course. However, our results should be insensitive to the inclusion 
of satellite quenching because satellites are typically too submm- 
faint to contribute to the (detectable) SMG population. 

Our method ensures that the mock sample has stellar mass 
functions, average specific SFRs, cosmic star formation history, 
and autocorrelation as a function of stellar mass that are consistent 
with observations. Because the model reproduces the mean SFR 
for a given stellar mass (including some scatter) and redshift rather 
than SFR(M*, z) for individual haloes, the elevation in SFR as- 
sociated with starbursts is not included, so we simply do not treat 
that subpopulation of SMGs in this work. This issue is discussed in 
detail in Sectionl4~5l 



2.2 Assigning submm flux densities 

To assign submm fluxes to the mock galaxies, we use the following 
fitting function, which was derived using the results of SUNRISE 
djonssonll2006l ; |jonsson et alfeOl O) dust radiative transfer calcula- 
tions performed on GADGET-3 dSpringelll2005l) smoothed-particle 
hydrodynamics simulations of isolated disc and merging galaxies 
(eq. 15ofH13): 



S 850 = 0.81 mJy 



SFR 



100 M Q yr- 



( M d 



V 10 8 M 



(1) 



where Ssso is the observed-frame 850- pm flux density and M<j is 
the dust mass (we describe the method for assigning Md below). 
This function can recover the submm flux of simulated galaxies in 
the redshift range z ~ 1 — 6 to within a scatter of 0.13 dex be- 
cause the negative K -correction makes the observed-frame Ssso of 



a fix ed galaxy SEP alm ost independent of z in the range z ~ 1 — 10 
(e.g jBlain et al.l2002l) . It under-predicts the flux at lower redshifts, 
significantly so for z < 0.5. However, because the normalization 
of th e SFR-Mt relation decreases strongly from z ~ 1 to z ~ 
(e.g jNoeske et alJl2007l ; lDaddi et al.ll2007l) and the volume in the 
range z = — 1 is a small fraction of the total z — — 6 vol- 
ume, both the ob served (e.g.,[Smolcic et al. 20121; lYun et al.ll20121 ; 



IWeiB et alJlioTl) and predicted dHavwardll2012l; H13) fraction of 
SMGs with 2 < 1 is small; thus, this unde r-prediction shoul d not 
significa ntly affect our results. Because the lEddingtonl d 1 9 1 31) bias 
(see also lHogg & Turnenll998l) can affect the shape of the number 
counts, we incorporate the scatter of 0.13 dex when assigning flux 
densities to the mock galaxies. 

The fitting function requires the SFR and M d of each mock 
galaxy . The SFR is already determined from the iBehroozi et alj 
( 120120 method|jTo assign dust mass values, we follow the method 
of HI 3, in which the gas fraction / a and metallicity Z are param- 
eterized as functions of stellar mass and redshift and it is assumed 
that 40 per cent of the metals are contained in dust grains. The func- 
tional form for / g (M*, z) is physically motivated, and the parame- 
ters ar e constrained to match z < 3 observations; see lHopkins et alj 
(2010) for details. The form for Z(M+, z) is based on the observed 
stellar mass-metallicity relation at different redshifts; see H13 for 
details. We note that because Sgso scales significantly sublinearly 



2 We have checked the effect of calculating the SFR following H 1 3 instead. 
For the majority of the sources, the SFRs and thus submm flux densities are 
very similar. Consequently, the results presented here are qualitatively in- 
sensitive to this choice, although there can be some quantitative discrepan- 
cies at the bright end because of the small number of sources that contribute 
there. 
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with both SFR and Md, the predictions are relatively insensitive 
to variations in these quantities; e.g., a factor of 2 change in either 
quantity changes Ssso by less than 50 per cent. 



2.3 Identifying blended sources 

Once submm flux densities are assigned to all mock galaxies, we 
perform a brute-force search for blended multiples. First, we se- 
lect only those sources with Sgso > 1 rnjy because fainter sources 
are below the flux density limits of current si ngle-dish surveys a nd 
interferometric follow-up observations (e.g.. iKarim et al.ll2012T) n 
This flux cut yields a minimum flux density of n mjy for n- 
component blended sources. Then, for each source i in this sub- 
sample, we compute the angular distance to all sources j in the 
subsample for which j > i (to avoid double-counting). If one or 
more sources is within a 'beam sizejj of source i, which we vary, 
we consider the sources to be blended and sum their Ssso values to 
calculate the total Ss50 f° r the blended source. 

For each blended source, we compute a measure of the redshift 
separation of the components by summing the redshift separations 
between the first component and all other components in quadra- 
ture: 



S S50 /2.3 = S U (mJy) 



1/2 



Az - J2^ 



Zl 



(2) 



where Zi denotes the redshift of the ith component and N is the 
total number of components of the blended source. This is a some- 
what arbitrary (yet natural) measure: for two components, Az is 
simply the difference in the redshifts of the two components. For 
sources with > 2 components, Az ~S> indicates that at least one 
of the components is at significantly different z than the others. Be- 
low, we define 'spatially associated' and 'spatially unassociated' 
blended sources by Az < 0.02 and Az ^ 0.02, respectively. The 
justification for and insensitivity of our results to the value of this 
cut are presented in Section [3~3] 

Fig.rjjpresents a simple illustration of the difference between 
spatially associated and unassociated projected pairs. The panels 
show the positions of a sample of mock galaxies projected onto the 
X — Y (left) and X — Z (right) planes. The galaxy pairs marked 
with red symbols are spatially associated in 3D (and thus also in 
projection). In contrast, the galaxy pairs marked with blue symbols 
appear to be associated in the X — Y projection, but the X — Z 



3 Sources fainter than the detection limit cause the zero-level of the submm 
map to be biased. However, this bias can be corrected w hen deriving num - 
ber counts from submm maps (see, e.g., section 3.2.3 of lWeifi et alj|2009ft . 
so we consider it reasonable to ignore the contribution of such sources to the 
predicted number counts. Using a different flux density limit in our model 
naturally changes the number of components of individual sources because 
for a lower (higher) flux density limit, the interferometric follow-up ob- 
servations would detect more (fewer) components. We have checked the 
effects of using alternate flux density limits of 5850 =0.5 and 2 mjy. The 
quantitative predictions for the absolute number counts change, but the rela- 
tive contributions of the subpopulations and the distributions of the redshift 
separations of components are relatively robust. Furthermore, the conclu- 
sions are qualitatively unchanged, and this uncertainty is subdominant to 
that associated with the treatment of blending. 

4 Throughout this work, we refer to the 'beam' and 'beam size' used in 
our model, but it should be understood that we are not claiming to precisely 
mimic the effects of blending (cf. fig. 1 of lHavward et alj|2012l) . Rather, 
our treatment is simply an approximation. This limitation is discussed in 
SectionR31 
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Figure 2. Cumulative number counts of model single-dish-detected SMGs 
(isolated discs and blends of multiple components) versus 5850 (mJy) for 
15-arcsec (red solid line) and 7-arcsec (blue dashed line) beams. The num- 
ber counts predicted when blending is not included are shown for compari- 
son (green dot-dashed line). All lines denote the median value taken over the 
eight mock catalogues, and the grey dashed error bars indicate the standard 
deviation. The black p o ints with solid erro r bars a r e observational data from 
Aretxaga et"aU feoill), lAustermann et al.1 feplOl). lHatsukade et al.l feOllI). 



Scott et all feOlOh. iKnudsen. van der Werf. & Kneitl <200g|) . IWeifi et alj 
(20091) , and lZemcov et al.N2010l) : see the legend for the corresponding sym- 
bols. The top axis indicates approximate 1.1-mm flux density values calcu- 
lated using Si.i ~ S85u/2.3 (see H13 for details). There are a sufficient 
number of blended SMGs to account for the observed Ssso Jj 12 mJy 
counts, which is in agreement with the recent suggestion that the brightest 
single-dish-detected SMGs are predominantly blended sources. 



projection clearly demonstrates that they are spatially unassoci- 
ated. Without accurate redshifts, an observer viewing these galax- 
ies along the Z axis cannot distinguish between these two types of 
projected pairs. This figure also illustrates how the spatial informa- 
tion provided by cosmological simulations is essential for treating 
the contribution of spatially unassociated projected multiples to the 
SMG population. 



3 RESULTS 

3.1 Cumulative number counts of blended and isolated-disc 
SMGs 

Fig.[2]shows the SMG cumulative number counts (deg~ 2 ) versus 
5*850 (mJy) predicted for beam sizes of 15 arcsec (which is rep- 
resentative of the resolution of the single-dish telescopes used to 
determine the observed counts) and 7 arcsec (for a conservative 
estimate of the effects of blending) and when blending is not in- 
cluded.The first two curves include blended SMGs and isolated disc 
SMGs that are not components of blended sources. All counts are 
the median counts taken over the eight mock catalogues, and the 
grey dashed error bars correspond to the standard deviation. The 
black points with solid error bars are observational data (see the 
caption for details). The simulation data plotted in this figure are 
presented in Table [T] Observed 1.1-mm counts have been approx- 
imately converted to 850-p.m counts using Ssso ~ 2. 35*1.1. This 
conversion is also used to show approximate Si.i values on the top 
axis; see H13 for details. 

Within the uncertainties, for the 15-arcsec beam, the pre- 
dicted counts of blended SMGs are consistent with most of the ob- 
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Table 1. Cumulative number counts for different amounts of blending. 



"S850" 


~ Si.i* 


iVl5<<(>S850) C 


N 7l ,(>s S50 y 


-Nnoblend(> SgSoY 


(mJy) 


(mJy) 


(deg- 2 ) 


(deg" 2 ) 


(deg" 2 ) 


2.0 


0.9 


3802 ± 206 


2362 ± 78 


731 ± 38.0 


3.0 


1.3 


2100 ± 185 


930 ± 64 


168 ± 12 


4.0 


1.7 


1093 ± 117 


379 ± 36 


51 ±4.7 


5.0 


2.2 


563 ± 76 


166 ± 23 


18 ±2.5 


6.0 


2.6 


288 ± 44 


80 ±15 


6.1 ±1.5 


7.0 


3.0 


151 ± 26 


38 ±7.5 


3.1 ±1.2 


8.0 


3.5 


85 ±15 


18 ±4.1 


1.5 ±1.0 


9.0 


3.9 


48 ±11 


12 ±2.2 


1.5 ±0.6 


10.0 


4.3 


29 ± 7.0 


6.1 ±1.9 


1.0 ±0.8 


11.0 


4.8 


17 ±4.9 


3.6 ±1.7 


5+ - 7 
U -°-0.5 


12.0 


5.2 


10 ±3.0 


2.6 ±1.0 


5+ - 6 
U '°-0.5 


13.0 


5.7 


5.6 ±2.4 


1.5 ±1.0 


< 0.5 


14.0 


6.1 


4.1 ±1.6 


1.0 ±0.6 


< 0.5 


15.0 


6.5 


2.0 ±0.8 


0.5 ±0.5 


< 0.5 



" 850-u.m flux density. Approximate 1.1-mm flux density calculated using Si.i ~ 
S850/2.3. c Cumulative 850-tim number counts for a beam size of 15 arcsec. d Same as 
column (3), but for a 7-arcsec beam. ' Cumulative number counts without the effects of 
blending included. 



served counts. The notable exception is the lAretxaga et al.l 1 1201 lh 
counts, which are sig nificantly higher than other measurements. 
lAretxaga et al.l ( 1201 lh argue that the origin of the discrepancy be- 
tween their counts and others observed for similarly sized fields 
is the impact of galaxy-galaxy weak lensing due to higher than 
average matter density along the line of sight. This hypothesis is 
supported by the association of the excess bright sources with fore- 
ground galaxy overdensities at 2 < 1.1. We do not include the 
effects of lensing here and thus would not produce this effect if it 
is the cause. However, the magnitude of the effect of lensing on the 
submm counts is still uncertain; see Section|4~5lfor further discus- 



example, for a ~ 1 deg 2 field, the flux density of the brightest non- 
blended source observed can vary in the range 5*850 ~ 10—12 mJy, 
and this effect is of course more significant for sma ller surveys. 
LESS, the survey from which the sources observed by Karim et al. 
were drawn, covers an area of 0.25 deg 2 , and the ISmolcic et alj 
sources were drawn from the central 0.7 deg 2 of the Cosmic Evo- 
lution Surve y (COSMOS ;IScoville et al.l2007r) field. Consequently, 
because the ISmolcic et al.l survey covered a larger area, the ex- 
pected flux density for th e brightest no n-blended SMG in their sur- 
vey is greater than for the lKarim et al.l survey. 



The difference between the predictions for the 15-arcsec beam 
(red solid line) and 7-arcsec beam (blue dashed line) indicates that 
the total number counts of the blended sources quantitatively de- 
pend on the beam size, as one naively expects. This suggests that 
the counts observed by telescop es of different reso lutions will dif- 
fer, as is already known (e.g.. ICasev et al.l 120131) . However, the 
difference also partially reflects uncertainty in the number counts 
predicted by our model because of our simple treatment of blend- 
ing (see Section |4.5t . Fortunately, there are multiple predictions 
that are robust to this aspect of the model. For example, for both 
beam sizes, the predicted counts agree within the uncertainties 
with the observed counts for 5sso > 12 mJy. Thus, the recent 
discover y that most SMGs with 5*850 ^ 12 mJy are blended 
sources JSmolcic et al.ll2012tlKarim et al.ll2012l) is consistent with 
the model presented here (but this result could be altered if star- 
bursts were included; see Section l4~T"l f or further discussion). 



Note also the significant uncertainty in the Sss 



9 mJy 



counts, which originates from sample variance and scatter in the 
SFR-halo mass relation. This affects all the counts, including the 
'no blending' case, because the maximum overdensity sampled 
varies with the field. It further affects the counts that include blend- 
ing by altering the probability that two high-sigma peaks have suf- 
ficiently small projected separation to be blended. Thus, even for 
independent fields with areas as large as ~ 2 deg 2 , the counts of the 
brighte st sources vary sign ifica ntly. This effec t could be the reason 
that the lKarim et alj <2012f) and lSmolcic et alj J2O120 samples vary 
in terms of the maximum flux density of non-blended SMGs. For 



3.2 Relative contributions of spatially associated and 
unassociated components 

Other robust predictions concern the relative contributions of spa- 
tially associated and unassociated components to the blended SMG 
population. Fig. [3] shows the fractional contribution of different 
subpopulations of blended SMGs to the total blended SMG cu- 
mulative number counts (i.e., the cumulative number counts of 
that subpopulation divided by the total cumulative number counts 
of blended sources) versus Ssso- The values are the medians of 
the individual ratios for the eight mock catalogues, and the grey 
dashed error bars indicate the standard deviations. We consider 
the spatially associateqj (solid lines) and spatially unassociated 
(dashed lines) subpopulations, which are defined by Az < 0.02 
and Az ^ 0.02, respectively. (The justification for this value is 
presented in Section 13.31 Here, we simply note that the results 
are insensitive to the specific cut because the Az distribution is 
strongly bimodal.) For both of the aforementioned subpopulations, 
the fractional contributions of blended sources composed of two 
(more than two) components are denoted with blue (red) lines, 



5 'Spatially associated' is a necessary but not sufficient condition for the 
individual components to be interacting galaxies. However, because a clas- 
sification based only on the difference(s) in the redshifts of the components 
will be most easily tested with forthcoming observations, we o pt to use 
this s imple criterion rather than a more sophisticated method (e.g.. lMorenol 
l20ii . 
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S 850 /2.3 ~ S, , (mJy) 



Table 2. Median Az values. 
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Figure 3. Relative contributions of various subpopulations to the cumu- 
lative number counts of blended SMGs for 15-arcsec (top) and 7-arcsec 
(bottom) beams. (Non-blended SMGs are ignored here.) The solid lines in- 
dicate blended sources for which all components are spatially associated 
(Az < 0.02), and the dashed lines are for blended SMGs in which at least 
one component is spatially unassociated (Az > 0.02). The blue (red) lines 
correspond to blended SMGs that are composed of 2 (> 2) components, and 
the (dashed) solid green line denotes the contribution of all (un)associated 
blended SMGs, regardless of the number of components. The lines indicate 
the median fraction at a given Sgso for the eight mock catalogues, and the 
grey dashed error bars indicate the standard deviation. At all flux densities 
plotted, spatially unassociated components contribute to ~ 30 — 80 per 
cent of the blended SMGs. For the 15-arcsec beam, the blended SMGs are 
dominated by > 2-component sources in which at least one component is 
spatially unassociated. The contribution of spatially unassociated compo- 
nents increases with the beam size. 



and the total fractional contributions of associated and unassoci- 
ated sources (i.e., independent of the number of components) are 
shown in green. 

One of the robust predictions of our model is that both spa- 
tially associated and unassociated blended sources contribute sig- 
nificantly at all flux densities and for the range of beam sizes ex- 
plored. Specifically, we expect that ~ 30 — 80 per cent of blended 
SMGs contain at least one component that is spatially unassoci- 
ated with the others. The fractional contribution varies with the 
flux density and beam size and may be sensitive to the treatment 
of blending, but the significance of the contribution is robust. For 
the 15-arcsec beam, the blended SMGs with Ssso > 3 mJy are 



Beam size 
(arcsec) Ssso > 3 mjy° Ssso > 7 mJy 6 



15 

7 



0.99 ±0.02 
0.88 ±0.04 



1.34 ±0.10 
1.11±0.17 



" Median and standard deviation of the Az values for 
each of the eight mock catalogues of blended SMGs 
with Ssso > 3 mJy and for the beam size specified in 
column (1). * Same as column (2), but for Ssso > 7 
mJy. 

dominated by > 2-component sources for which at least one com- 
ponent is spatially unassociated with the others. In contrast, for the 
7-arcsec beam, the typical blended SMG in our model is composed 
of two associated galaxies, but at all flux densities shown, at least 
~ 30 per cent of the blended SMGs have one or more spatially 
unassociated components. 

Fig. [3] also indicates that the relative contribution of spatially 
unassociated sources is greater for the 15-arcsec beam. This is sim- 
ply because the greater beam size increases the likelihood of mul- 
tiple physically and spatially unassociated galaxies falling with the 
same beam. In contrast, the relative contribution of spatially associ- 
ated blended SMGs is less for the 15-arcsec beam for two reasons: 
1. The number of additional physically associated galaxies blended 
when the beam size is 15 arcsec is relatively small because the pro- 
jected separation of such systems is typically less than the 7-arcsec 
beam (i.e., ~ 65 kpc at z ~ 2 — 3); thus, the 7-arcsec beam is suf- 
ficient to blend these sources. 2. Some sources that are spatially as- 
sociated for the 7-arcsec beam become classified as spatially unas- 
sociated for the 15-arcsec beam because one or more additional, 
spatially unassociated galaxies become blended with the spatially 
associated sources when the beam size is increased to 15 arcsec. 
For the same reason, the number of sources with > 2 components 
is significantly greater when the beam size is 15 arcsec. 

Finally, Fig. [3] indicates that for the 7-arcsec beam, the frac- 
tional contribution of spatially associated blended sources increases 
with flux density. We speculate that the cause of this trend is the 
following: to have the necessary Ssso, a galaxy must be very mas- 
sive and thus correspond to relatively rare, high-overdensity peaks. 
The required mass increases with flux density. Thus, the higher- 
5850 sources trace increasingly rare peaks, and the likelihood that 
two or more sufficiently rare peaks have projected separation less 
than 7 arcsec decreases with Ssso- Consequently, t he brightest 
blende d sou rces correspond to protoclusters (see, e.g. jDaddi et alj 
l2009al lbl and llvison et alj|2013l for real-Universe examples of such 
SMGs) rather than chance projections of multiple extremely mas- 
sive galaxies because massive haloes clu ster more strongly, and 
this effect increases with red shift (e.g., ISheth & Tormenl 11993 : 
lAngulo. Baugh. & Lacevll2008l) . 



3.3 Redshift separations of the components of blended SMGs 

The final testable prediction that we make concerns the distribution 
of the redshift separation of the individual components of blended 
SMGs. Fig.|4]shows the distributions of log Az for blended sources 
that are brighter than two different flux cuts, Ssso > 3 (solid lines) 
and 7 (dotted lines) mJy for both beam sizes (15 arcsec in red and 7 
arcsec in blue). The distributions are strongly bimodal (we interpret 
this bimodality below), and there is a clear minimum at Az ~ 0.02 
(marked by the vertical dashed line), which is why we use this 
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Figure 4. The solid (dotted) lines show the distribution of the logarithm 
of the redshift separation of the individual components (as defined in Eq. 
[2} of the blended submm sources in one of the mock lightcones with total 
blended flux density Ssso > 3 ( 7 ) mJy for both the 15-arcsec (red) and 
7-arcsec (blue) beams. The distribution is strongly bimodal, which justi- 
fies the Az cut used to separate the spatially associated and unassociated 
subpopulations; this cut (Az = 0.02) is indicated by the vertical dashed 
line. The median values are Az ~ 0.9 — 1.3, depending on the beam size 
and flux cut (see Table[2). Note that the shapes of the distributions are rela- 
tively insensitive to the beam size, although for the 7-arcsec beam, a higher 
fraction of the sources are spatially associated. 



value as the cut to distinguish the spatially associated and unas- 
sociated subpopulations. (Using a smaller, more physically mean- 
ingful cut, such as Az — 0.01, does not significantly alter the re- 
sults because of the strong bimodality.) As discussed above, for the 
7-arcsec beam, a higher fraction of the sources are spatially associ- 
ated. 

Many of the galaxies in the spatially associated subpopulation 
(Az < 0.02) are also physically associated, i.e., gravitationally 
bound galaxies that will merge. The galaxy-pair SMG s (early-stage 
mergers) discussed in detail in lHavward et al.l ( 120120 and H13 are 
a subset of the spatially associated galaxies. This category also in- 
volves galaxies that are undergoing fly-by encounters but will not 
merge. Typically, the spatially associated galaxies are in distinct 
haloes, and the separations can be as large as ~ 10 Mpc (comov- 
ing). The spatially unassociated galaxies (Az ^ 0.02) include at 
least one component that is simply a chance projection, i.e., is com- 
pletely unrelated with the other(s). The separations amongst the 
components of the spatially unassociated galaxies range from a few 
Mpc to the box size. 

The median and standard deviation of the Az values (calcu- 
lated from the individual values for each mock catalogue) for the 
spatially unassociated subpopulations are listed in Table [2] Note 
that the median values are of order 1, which is a natural conse- 
quence of the redshift distribution of SMGs in our model: most 
SMGs are in the range z ~ 2 — 4, and the median redshift is ~ 3[j 
Thus, the typical separation of components is naturally ~ 1. The 



median value of Az is relatively insensitive to the beam size, flux 
density cut, and mock catalogue used; thus, it is a robust, testable 
prediction of our model. The maximum Az varies significantly for 
the different mock catalogues and ranges from z ~ 4 — 6. 



4 DISCUSSION 

4.1 Can the brightest SMGs be explained via blending alone? 

The recent obse rvational evidence ( IKarim et alj 120121 : 
ISmolcic et"al] 1201 2h that most of the single-dish SMGs with 
5*870 > 12 mJy are blends of multiple components was part of 
the impetus for this work. In particular, we sought to address 
whether our model could predict a sufficient number of blended 
SMGs to account for the observed number counts of sources 
with 5870 ^ 12 mJy because in the H13 model, such sources are 
predominantly starbursts (the models are compared in detail in the 
next subsection). 

By including the contribution from spatially unassociated 
components (i.e., chance projections), we predict blended SMG 
number counts (for a 15-arcsec beam) that are consistent with 
the observed count s within the uncertainties (except for the 
lAretxaga et alj|201 ll counts, which are thought to be significantly 
boosted by weak lensing, as discussed above). Even for the more 
conservative treatment of blending, in which we use a 7-arcsec 
beam, the predicted and observed counts are not significantly dif- 
ferent for SWo > 12 mJy. Consequently, assuming that our the- 
oretically derived scaling relation is approximately accurate, there 
are a sufficient number of blended SMGs to explain the 5*870 ^ 12 
mJy counts. 

However, because the model presented here does not include 
merger-induced starbursts (see Section 14.5) , we cannot yet con- 
clusively argue that the 5s7o > 12 mJy SMGs are exclusively 
blends of multiple components. Furthermore, we must account for 
the sources with Ss70 > 12 mJy that a re n ot blended sources, be- 
cause some examples are known (see iHavwardl 120131 for further 
discussion). Both issues are deferred to future work, which we out- 
line below. 



4.2 The importance of mock observations for testing models 

We have seen that a given single-dish-detected SMG can actually 
correspond to multiple physical galaxies|j However, it is often as- 
sumed that an SMG corresponds to a single (real or model) galaxy. 
Although this assumption is valid for isolated discs and reasonable 
for late-stage mergers, it is questionable for early-stage mergers, in 
which the progenitor discs are well separated and not yet strongly 
interacting, and definitely invalid for projected multiples. This er- 
ror, which is commonly committed by theoreticians but better un- 
derstood amongst observers, can be partially attributed to the term 
'submillimetre galaxy', which is somewhat misleading|j 

Were this simply a matter of terminology, it would not be 



6 At z < 2, the SFRs and (thu s submm fluxes) of star-forming galax- 
ies are significantly lower (e.g., [ Noeske et al.l. 2 007; Daddi et al. 2007; 
iMagnelli et"aI]|2009Ll20ilT;lKarim et al.ll201ll ; IWhitaker et alJboTa, a s 



nificant frac ti on of massive gal a xies have been quenched (e.g., Befle^al 



2004 [2007; Faber et al. 200^ lllbert et al.1 l2010l 120131; iDavidzon et ~ 



20131) , and the volume probed is relatively small; all of these reasons cause 



the SMG counts to drop significantly at z < 2. At z > 4, the scarcity of 



massive galaxie s and decreased dust content of galaxies causes the counts 
to decrease. See Hay ward 2012 and H13 for more details. 

7 By 'galax y', we refer to a v i rialize d dark matter halo that contains 
baryons; see IWillman & S trader (2012|) for a more observationally moti- 
vated definition. 

8 Perhaps 'submillimetre source' is a better term, but the acronym 'SMS' 
is already widely used for another purpose. 
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a concern, because the astronomical literature is rife with mis- 
leading terminology. However, equating SMGs with individual 
'galaxies' in simulations or semi-analytical models (SAMs) can 
lead to qualitatively inaccurate conclusions because of the signif- 
icant differences between predictions for individual model galax- 
ies (e.g., the 'no blending' curve in Fig. [2]l and predictions for 
model SMGs (e.g., the other curves in Fig.|2j; to make predictions 
that are classified into the latter category, a model must attempt 
to at least crudely mimic observations. This model-c omparison 
problem is n ot specific to SMGs , of course (see, e.g., Lotz et alj 
20081 l2010djbl; IWuvts et alj|2009l. l20ld : IScannapieco et alj|20ld : 
Snyder et alj201 Ul2012l) . but it is especially acute for them because 
1. observations of SMGs are more limited than for local and/or less 
dusty galaxies and 2. the small dynamic range of submm flux den- 
sity currently probed and the steepness of the number counts im- 
plies that relatively small changes in flux density - caused by, e.g., 
blending - can qualitatively affect the conclusions. 



4.3 Implications for the inferred SFRs of SMGs 

Properly treating blending may reconcile the longstanding discrep- 
ancy between the far-IR and (sub)mm counts predicted by cosmo- 
logical simulations and SAMs and those observed. Furthermore, 
the ubiquity of blended submm sources has an important impli- 
cation for inteipreting SMG observations: because many of the 
brightest sources are actually multiple unrelated galaxies, the SFRs 
inferred under the assumption that a single submm source corre- 
sponds to a single galaxy overestimate the true SFRs of individual 
galaxies. This overestimate of the abundance of the most rapidly 
star-forming galaxies may relieve the tension betwee n the SFRs 
of the mock SMGs in the cosmological simulation of iDave et alj 
1 120 101) and those inferred from single-dish observations. 

To obtain the true SFRs of individual galaxies, resolved pho- 
tometry of the individual components is required. In the UV- 
optical, this process is limited by the significant attenuation that 
is characteristic of SMGs. In the IR-(sub)mm, the resolution is the 
limitation. The A LMA continuum m easurements from the LESS 
follow-up survey JKarim et alj|2012c Hodge et al., submitted) will 
help in this regard, but additional high-resolution data in the far- 
IR, particularly shortward of the peak of the dust emission, are 
necessary to accurately constrain the IR luminosity and thus infer 
the SFRs of the individual components. Therefore, the question of 
whether the SFRs of galaxies in cosmological simulations and real 
galaxies are in tension should be revisited when possible. 



4.4 Comparison with the H13 model 

Although the model presented in this work and that of H13 have 
some ingredients in common, they differ in the following impor- 
tant aspects: (1) Here, we directly include information about large- 
scale structure by using mock galaxy catalogues generated from the 
Bolshoi simulation. In HI 3, the cosmological context was obtained 
solely through the use of merger rates calibrated to cosmological 
simulations and an assumed stellar mass function. Consequently, 
no spatial information was included, and thus chance projections 
of spatially unassociated galaxies could not be treated. (2) In HI 3, 
the effects of blending were included for early-stage mergers, but 
in a cruder manner (by assuming a fixed timescale) than in this 
work. (3) Because the procedure to generate the mock galaxy cat- 
alogues relies on an SFR-halo mass relation (that depends on red- 
shift and includes some scatter), the full SFR enhancement caused 



by merger-induced starbursts is not included in this work. In HI 3, 
the submm light curves for merger-induced starbursts were directly 
calculated by performing dust radiative transfer on hydrodynamical 
simulations. Consequently, the SFR enhancement in such objects 
was treated in the H13 model. 

The differences described above imply that neither the H13 
model nor that presented here can fully treat all proposed SMG 
subpopulations. Thus, the predictions of the two models cannot 
be directly compared. However, additional insight can be obtained 
by combining some of their results. Here, we have demonstrated 
that blending of unassociated galaxies can significantly boost the 
submm number counts. Thus, our results reinforce the conclu- 
sion of H13 that given a demographicalfy accurate model galaxy 
population (i.e., one with stellar mass functions and merger rates 
consistent with observations), the observed submm counts can be 
matched without recourse to variation in the stellar initial mass 
function. Furthermore, if spatially unassociated multiples were in- 
cluded in the H13 model, it is possible that the model would over- 
predict the observed counts. Such an overprediction could imply 
that the normalisation of the stellar mass function used in the H13 
model is too high or that the H13 simulations are missing some 
important form of feedback that limits the SFR in starbursts. 

4.5 Limitations of our model 

The strength of the model presented in this work is that it naturally 
includes information about large-scale structure and thus enables 
us to physically characterise the effects of blending for both spa- 
tially associated and unassociated galaxies; in contrast, H13 treated 
only the blending of the components of early-stage mergers, and 
the treatment was significantly cruder. However, the method used 
here naturally also has some limitations. 

First, the method by which we assign SFRs to galaxies treats 
scatter in the SFR-halo mass relation at a given mass and redshift, 
but it does not faithfully reproduce the SFR evolution of individual 
haloes. In particular, the elevation in SFR associated with merger- 
induced starbursts (i.e., the high-SFR tail at a given halo mass), 
which is genet ically predicted by idealised simulations of major 



mer gers (e.g.. | Hernquistjll989t iBarnes & Hernquistlll99U 



Mihos & Hemquist 1996; Springel, Di Matteo, & Hernquist 



1996; 



2005; 



ICox et alj|2006l ; I Younger et al.ll2009h . is not included. It would be 
possi ble to include this effec t in an empirically motivated manner 
(e.g.. iBethermin et alj 120121) . However, because merger-induced 
starbursts were treated extensively in HI 3, and the primary focus of 
this work is to investigate the contribution of spatially unassociated 
components to blended SMGs, we opt to simply neglect merger- 
induced starbursts here. 

Subsets of the spatially associated 2-component and both spa- 
tially associated and unassociated > 2-component blended SMGs 
should be undergoing merger-induced starbursts. (The reason that 
some of the spatially unassociated > 2-component blended SMGs 
should be undergoing starbursts is that in those sources, a subset of 
the components can be spatially and physically related.) Because 
we do not include the elevation in SFR from starbursts, we under- 
estimate the submm fluxes of such objects. However, because the 
duty cycle of the starburst phase is significantly less than that of the 
'galaxy-pair SMG' phase (H13) and not all encounters cause star- 
bursts, the fraction of blended SMGs that are undergoing starbursts 
should be relatively small. For more typical spatially associated 
multiples, which have large separations (i.e., > 10 kpc), the SFR 
elevation should be relatively modest (e.g., IScudder et alj 120121 : 
lLanz et alj2013D . Consequently, our conclusions regarding the rel- 
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ative contributions of spatially associated and unassociated compo- 
nents to blended SMGs and the distributions of the redshift sepa- 
rations of the components would not be significantly affected by 
the inclusion of starbursts (but the quantitative predictions regard- 
ing, e.g., the number counts should be affected, and the conclusions 
regarding the brightest sources may differ). Still, there is strong ob- 
servational evidence that galaxies above the main sequence (i.e., 
merger-induced starbursts and gravitationall y unstable disc galax- 
ies) contribute to the SM G population (e.g.. lMagnelli et alj|2012l ; 
iMichalowski et alj|2012h . and H13 predicted a significant starburst 
contribution to the SMG population. Thus, we should and will ad- 
dress this shortcoming of the model presented here in future work. 

Furthermore, the treatment of blending is still somewhat 
crude, albeit significantly less crude than that used in H13 and far 
superior to simply ignoring blending. The difference between the 
predicted counts for the 15- and 7-arcsec beam sizes illustrated by 
Fig.|2]indicates that the treatment of blending can significantly af- 
fect some predictions of the model. Although there should certainly 
be differences in the blended SMG populations identified using 
telescopes of different beam sizes, the predictions for a given beam 
size in the model should not be considered to exactly correspond to 
observations from a telescope of the same beam size. Rather, for the 
typical single-dish SMG surveys, which have been performed with 
telescopes with ~ 15-arcsec beams, our predictions for the 15- and 
7-arcsec beams should be considered somewhat liberal and very 
conservative estimates of the number counts of blended SMGs, re- 
spectively. 

Finally, we have neglected the effects of gravitational lens- 
ing on the predicted n umber counts. As already mentioned above, 
I Aretxaga et al.l ( 120 1 ll) suggest that galaxy-galaxy lensing is the rea- 
son that their counts are significantly higher than those for other 
fields (see also lAustermann et alj2009i) . Furthe rmore, a populatio n 
of strongly lensed SMGs has been observed dVieira et al.l 1201 01) . 
However, the effect of lensing for sources with 5*850 % 10 mjy is 
poorly constrained: estimates of the typic al magnification for such 
sources ra nge from ~ 10 per cent (e.g., [Pa ciga, Sco tt, & Chapinl 
l2009t iLima et al.ll2010bl ; lLima. Jain, & Devlirj|2010al) to of order 
10 JHarris et alj|2012r) . Thus, we have opted to defer a treatment of 
lensing to future work. 



4.6 Future work 

In future work, we will combine our fitting functio n for submm flux 
densit y wi th an updated version o f the SAM of ISomerville et alj 
( 120081) and lSomerville et al.l 020120 to produce mock submm maps 
of the galaxies from the SAM. Then, to make predictions for var- 
ious telescopes, we will convolve the maps with the appropriate 
beams. This approach will enable us to simultaneously treat the iso- 
lated disc, merger-induced starburst, and both spatially associated 
and unassociated blended SMG subpopulations in order to present 
a more complete picture of the SMG population. Furthermore, we 
will investigate the effects of modifying the prescriptions in the 
SAM for, e.g., star formation and feedback on the mock SMG pop- 
ulation. Finally, the spatial information available in our mock SMG 
catalogues makes them a natural tool for investigating SMG clus- 
tering and the effects of lensing on submm number counts. 



5 CONCLUSIONS 

We created mock SMG catalogues by generating lightcones from 
the Bolshoi simulation, assigning SFRs and dust masses to the 



mock galaxies in an empirically constrained manner, and using a 
fitting function from previous work to assign submm flux densities. 
Then, the mock SMG catalogues were used to investigate the ef- 
fects of blending of multiple galaxies on the SMG population. Our 
principal conclusions are the following: 

(i) For 5*850 > 12 mjy, the blended SMG number counts pre- 
dicted by our model are similar to those observed, which is consis- 
tent with the recent observational evidence that most of the bright- 
est SMGs are blended sources. 

(ii) The d ifference in the maxim um flux density of non -blended 
SMGs in the lSmolcic et alj J2012I) and lKarim et alj 1201 2h samples 
is consistent with the variation in the flux density of the brightest 
non-blended SMG amongst our mock SMG catalogues. Thus, the 
difference may simply be a consequence of sample variance and 
scatter in the SFR-halo mass relation. 

(iii) Both spatially associated and unassociated blended SMGs 
contribute significantly to the population: ~ 30 — 80 percent of the 
blended SMGs in our model contain one or more components that 
are spatially distinct from the other(s). The fractional contribution 
depends on both the flux density and beam size. For a 15-arcsec 
beam, blends of > 2 galaxies in which at least one galaxy is spa- 
tially unassociated with the others dominate the blended sources 
with total 5850 > 3 mJy. 

(iv) The distribution of the redshift separations of the compo- 
nents of blended SMGs is strongly bimodal. For the spatially unas- 
sociated subpopulation, the typical redshift separation is Az ~ 
0.9 — 1.3, depending on the beam size and flux cut. 

We stress that the last two conclusions are bona fide predic- 
tions in the sense that (to our knowledge) currently available obser- 
vational data are insufficient to test them, but they will be tested as 
soon as spectroscopic or sufficiently accurate photometric redshifts 
of the individual components of a sufficient number of blended 
SMGs are available. 
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